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G^ ■ ABSTRACT 

O . 

OO , Aims. The BL Lac object S5 0716+71 was observed in a global multi-frequency campaign to search for rapid and correlated flux den- 

. sity variability and signatures of an inverse-Compton (IC) catastrophe during the states of extreme apparent brightness temperatures. 

Methods. The observing campaign involved simultaneous ground-based monitoring at radio to IR/optical wavelengths and was cen- 
tered around a 500-ks pointing with the INTEGRAL satellite (November 10-17, 2003). Here, we present the combined analysis and 
Jv> , results of the radio observations, covering the cm- to sub-mm bands. This facilitates a detailed study of the variability characteristics 

j_j ■ of an inter- to intra-day variable IDV source from cm- to the short mm-bands. We further aim to constrain the variability brightness 

' temperatures (Tg) and Doppler factors (6) comparing the radio-bands with the hard X-ray emission, as seen by INTEGRAL at 3- 

200 keV. 

Results. 0716+714 was in an exceptionally high state and different (slower) phase of short-term variability, when compared to the 
past, most likely due to a pronounced outburst shortly before the campaign. The flux density variability in the cm- to mm-bands is 
dominated by a ~ 4 day time scale amplitude increase of up to ~ 35 %, systematically more pronounced towards shorter wavelengths. 
The cross-correlation analysis reveals systematic time-lags with the higher frequencies varying earlier, similar to canonical variability 
on longer time-scales. The increase of the variability amplitudes with frequency contradicts expectations from standard interstellar 
scintillation (ISS) and suggests a source-intrinsic origin for the observed inter-day variability. We find an inverted synchrotron spec- 
trum peaking near 90 GHz, with the peak flux increasing during the first 4 days. The lower limits to Tb derived from the inter-day 
variations exceed the 10'^ K IC-limit by up to 3-4 orders of magnitude. Assuming relativistic boosting, our different estimates of 
6 yield robust and self-consistent lower limits of 5 > 5 - 33 - in good agreement with Svlbi obtained from VLBI studies and the 
IC-Doppler factors S/c > 14-16 obtained from the INTEGRAL data. 

Conclusions. The non-detection of S5 0716+714 with INTEGRAL in this campaign excludes an excessively high X-ray flux associ- 
ated with a simultaneous IC catastrophe. Since a strong contribution from ISS can be excluded, we conclude that relativistic Doppler 
boosting naturally explains the apparent violation of the theoretical limits. All derived Doppler factors are internally consistent, agree 
with the results from different observations and can be explained within the framework of standard synchrotron- self-Compton (SSC) 
jet models of AGN. 
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1. Introduction 

Rapid intensity and polarization variability on time scales of 
hours to days is frequently observed in compact blazar cores 
at centimeter wavelengths. Since its discovery in 1985 such 
IntraDay variability (IDV, Witzei et al. 1986, Heeschen et al. 
1987) has been found in a significant fraction (~ 20-30%) of 
flat-spectrum quasars and BLLac objects (Quirrenbach et al. 
1992, Kedziora-Chudczer et al. 20 01, Lovell e t al. 2003). The 
'classical' IDV sources (of type II. IHeeschen e t al., 1987) show 
flux density variations in the radio bands on time scales of 
< 0.5 -2 days and variability amplitudes ranging from a few up to 
30%. Stronger an d often faster variability i s seen in both linear 
polarization (e.g. lOuirrenbach et al.l 119891: iKraus et aT, 20031) 
and m ore recently also in circular polarization (Macquart et aU 
I2OOOI) . 

When the observed IDV is interpreted as being source- 
intrinsic, the apparent variability brightness temperatures of 
Tg > 10 '^ K, largely exceed the inv erse Compton (IC ) limit 
of 10'^ K jKellermann & Pauhnv-Toth, 1969: ReadheaA 11994 ) 
by se veral orders of magnitude (see also IWagner & Witzei , 
Il995h . Consequently, extreme relativistic Doppler-boosting fac- 
tors {6 > 100) would be required to bring these values down 
to the IC-limit. As such high Doppler-factors (and related bulk 
Lorentz factors) are not observed in AGN jets with VLBI, al- 
ternate jet (and shock-in-jet) models use non-spherical geome- 
tries, which allow for additional r elativistic corrections (e.g. 
ISpada et al.1 119991: loian et all 1199 lb . It is also possible that the 
brightness temperatures are i ntrinsically high and coherent emis- 
sion processes are involved dBenfordl [19921 iBenford & Leschl 
Il998l) . Another possibility to explain intrinsic brightness tem- 
peratures of > 10'^ K is a more homogeneous ordering of the 
magnetic field (lOian et al.l l2006). 

Owing to the small source sizes in compact flat-spectrum 
AGN, interstellar scintillation (ISS) is an unavoidable process. 
For a small number of more rapid, intra-hour variable sources, 
seasonal cycles and/or variability pattern arrival time delays are 
observed. This is used as convincing argument that such rapid 
variabihty is caused by ISS (PKS 0405-385, PKS 1257-326 
and J1819-H3845. e.g. iJauncev et all [2000l iBignall et all [200I 
[Dennett- Thorpe & de BruvnT 20021) . The situation for the slower 
and 'classical' (type II) IDV sources is less clear The much 
longer and complex variability time scales and the existence 
of more than one characteristic variability time-scale, did not 
yet allow to unambiguously establish seasonal cycles nor ar- 
rival time measurements and by one of this proof that IDV of 
type II is solely due to ISS. In fact, many of these sources 
show an overall more complex variability behavior, (e.g. QSO 
0917-H624, BL 0954-H658, J 1 128-H5925: iRickett et aLl [2001 ; 
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Jauncev & Macquart, 2001; Fuhrmann et al., 2002^, IFuhrmann , 
2004; Bernhart et al., 2006; Gabanyi et al., 2007), which points 
towards a more complicated and time-variable blending between 
propagation induced source extrinsic effects and source intrinsic 
variability. 

In order to shed more light on the physical origin of IDV 
in such type II sources, we performed a coordinated broad- 
band variability campaign on one of the most prominent IDV 
sources. Since the late 1980's, the compact blazar S5 0716+71 
(hereafter 0716+714) showed strong and fast IDV of type II 
whenever it was observed. In the optical band, 0716+714 is 
identified as a BLLac, but has unknown redshift. A redshift 
of z >0.3 was deduced using the limits on the surface bright- 
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ll993l lWagner et al.lll996tlSbarufatti et an.l2005l [Nilsson et al. 



20081). In this paper we will use z >0.3, which is within 
the measurement uncertainty of the previous redshift estimates. 
Consequently, only lower limits to any linear size and bright- 
ness temperature measurement can be obtained. The observed 
IDV time scales imply - via the light travel time argument 
- source sizes of only light-hours (corresponding to micro- 
arcsecond scales) leading to lower limits of Tb > 10'^ K. Direct 
measurement of the nucleus with space- VLBI at 5 GHz gives 
a robust lower limit of Tb > 2 x 10'^ K, implying a min imum 
equipartition Doppler factor of <5 > 4 (B ach et aU l2006l) . The 
VLBI kinematics leads to estimates of the bulk jet Lorentz factor 
Tmin ^ 16-21, based on the measured jet speeds dJorstad et al.1 
I2OOII [Bach et al.1 12005|) . This is is unusually high for a BL Lac 
object. 

0716+714 is one of the best stud i ed blazars in the sky (e.g . 
Wagner et al.1 1 19961 [pabuzda et al.1 ll998l lOuirrenbach et all 
2OOOI iKraus et all l2003l iBach et all I200I |2006|) . It has been 
observed repeatedly in various multi-frequency campaigns 
(e,g^ Wagner etal., 1990, 1996; Otterbein & Wagner, 199^ 
iGiommi et aUll999llTaghaferri et al.ll2003llRaiteri et al.ll2003[) 
and is well known to be extr emely variable (< hours to months) 
at radio to X-ray bands (e.g. IWagner et al.l 1 19961 iRaiteri et al 



12003 , and references therein). 0716+714 is so far the only source 
in which a correl ation between the radio/optical variability 
was o bserved (see lOuirrenbach et all Il99ll IWagner & Witzell 
1 19951) indicating that the radio/optical emission has a com- 
mon and source-intrinsic origin, during the observations in 1990 
(lOian et all [19961120021) . Further, the detection of IDV at mm- 
wavelengths and the obse rved frequency dependence of the 
IDV variability amplitudes dKrichbaum et al.ll2002tlKraus et all 
2003) make it difficult to interpret the IDV in this source solely 
by standard ISS. 

In a source-intrinsic interpretation of IDV it is unclear if 
and for how long the IC-limit can be violated dSlvshl \l99% 
lKellermannll2002l) . The quasi-periodic and persistent violation 
of this limit (on time scales of hours to days) would imply 
subsequent Compton catastrophes, each leading to outbursts of 
IC scattered radiation. The efficient IC-cooling would rather 
quickly restore the local Tb in the source, and the scattered radi- 
ation s hould be observable as enhanced em i ssion in the X-/y -ray 
bands (Kellermann & Pauli nv-TotR 1 19691 iBloom & Marschen, 
1996; Ostorero et al., 2006). In order to search for multi- 
frequency signatures of such short-term IC-flashes, 0716+714 
was the target of a global multi-frequency campaign carried out 
in November 2003, and centered around a 500-ks INTEGRAlE 
pointing on November 10-17, 2003 ('core-campaign'). In or- 
der to detect contemporaneous IC-limit violations at lower ener- 
gies, quasi-simultaneous and densely time-sampled flux density, 
polarization and VLBI monitoring observations were organized, 
covering the radio, millimeter, sub-millimeter, IR and optica^ 
wavelengths. 

Early results of this campaign including radio data at two 
frequencie s, optical and I NTEGRAL soft y-ray data were pre- 
sented by lOstorero et al.l ^)06), who also showed the simul- 
taneous spectral energy distribution. The results from the mm- 
observations (3 mm, 1.3 mm) perf ormed with the IR AM 30 m 
radio-telescope were presented bv lA gudo et al.l (1200 6). The re- 
sults of the VLBI observations will be presented in a forthcom- 
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ing paper. In this paper III we present the intensity and polar- 
ization data obtained with the EfFelsberg 100 m radio telescope 
(5, 10.5, & 32 GHz) during the time of the core-campaign. We 
then combine these data with the (sub-) millimeter and optical 
data in order to determine the broad-band variability and spectral 
characteristics of 0716-1-714, with special regard to (i) the intra- 
to inter-day cm-/mm-variability, (ii) the brightness temperatures 
and IC-limit, and (iii) the Doppler factors combining the radio 
and high energy INTEGRAL data from this campaign. 

2. Observations and Data Reduction 

2.1. Participating Radio Observatories 

The intensive flux density monitoring of 0716+714 was carried 
out at the time of the core campaign between November 08-19, 
2003 (JD 2452951.9-2452962.8). In total, 7 radio observatories 
were involved, covering a frequency range from 1 .4 to 666 GHz 
(wavelengths ranging from 21 cm to 0.45 mm). In Table [l] the 
participating observatories, their wavelength/frequency cover- 
age and observing dates are summarized. 

At the Effelsberg telescope, 0716+714 and a number of sec- 
ondary calibrators were observed continuously at 5, 10.5, & 
32 GHz, with a dense time sampling of about two flux den- 
sity measurements per hour, source and frequency. The spe- 
cial ID V observing strategy and matched data reduction pro- 
cedure dKraus et al I l2003h enabled high precision flux density 
and polarization measurements and facilitate the study of rapid 
variability on time scales from about 0.5 hrs to 7 days. At 3.5 
& 1.3 mm wavelengths with the IRAM 30 m-telescope on Pico 
Ve leta, a similar observing strategy was applied and is described 
bv lAgudo et al.l ( l2006l) . At Effelsberg polarization was recorded 
at ^6 cm (4.85 GHz) and ^2.8 cm (10.45 GHz), at the IRAM 
30 m-telescope at 3.5 mm (86 GHz). The matched observing and 
calibration strategies at Effelsberg and Pico Veleta now allow for 
the first time a detailed study (and e.g. cross-correlation analysis) 
of the short-term variability of an IDV source from centimeter to 
short-millimeter wavelengths. 

Complementary, but less dense in time sampled flux density 
measurements were also performed with the telescopes at Kitt 
Peak (12 m), Mt. Graham (SMT/HHT 10m), and Mauna-Kea 
(JCMT, 15m), which covered the higher frequency bands (86 to 
666 GHz, see Table 1). The Westerbork interferometer provided 
a few single measurements, which extends the frequency cover- 
age to 1 .4 and 2.3 GHz. The light curves obtained at 32 GHz with 
the 100 m telescope and at 37 GHz wit h the Metsahov i 14 m - 
telescope have been already shown by lOstorero et alj (l2006h . 
The analysis of the 32 GHz Effelsberg data has been improved. 



ments using cross-scans (azimuth/elevation direction) with the 
number of sub-scans matching the source brightness at the given 
frequency ( see Heeschen et al., 1987; Ouirrenbach et al, 199^ 
iKraus et all [20031 for details). Depending on the number of 
sub-scans (typically 4, 8 or 12), a single cross-scan lasted 2- 
4 minutes. A duty cycle consisted of subsequent observation 
of the target source and several nearby secondary calibrators 
(steep-spectrum, point-like, non-IDV sources, e.g. 0836+710 
and 0951+699). Within such a duty-cycle, each source was ob- 
served consecutively at all three frequencies (in the order 32, 
10.5 and 4.85 GHz). During the 7 days observing period and 
our continuous 24hr/day time coverage, we obtained for each 
source about 2 measurements per hour and frequency. Primary 
flux density and polarization angle calibrators (e.g. 3C 286, see 
Table O were observed every few hours. Among our target 
source 0716+714, two additional IDV sources were observed: 
0602+673 and 0917+62^ 

The data reducti on was done in the standard m a nner, and 
is described in e.g. iKraus et aL 1 (120031) : iFuhrmannI ( |2004|) . It 
consists of the following steps: (i) baseline subtraction and fit- 
ting of a Gaussian profile to each individual sub-scan; (ii) flag- 
ging of discrepant or otherwise bad sub-scans/scans; (iii) cor- 
recting the measured amplitudes for small residual pointing er- 
rors of the telescope (< 5-20 arcseconds); (iv) sub-scan averag- 
ing over both slewing directions (the Q/U data and the A 9 mm 
total power data were averaged before Gaussian fitting to en- 
hanced the signal-to-noise ratio and the quality of the Gaussian 
fits); (v) opacity correction using the Tsyj-measurements ob- 
tained with each each cross-scan; and (vi) correction for re- 
maining systematic gain-elevation and time-dependent gain ef- 
fects using gain-transfer functions derived from the dense in 
time-sampled secondary and primary calibrator measurements 
assuming their stationarity. Finally, the measured antenna tem- 
peratures for each observed source were linked to the absolute 
flux-density scale, using the frequent p rimary calibrator mea- 
surements (Baars eTall Il977l; lOtt et all IT994) . Their flux den- 
sities and polarization properties are summarized in Table[3] 

The individual flux density errors are composed of the statis- 
tical errors from the data reduction process (including the errors 
from the Gaussian fit, the weighted average over the sub-scans, 
and the gain and time-dependent corrections) and a contribution 
from the residual scatter seen in the primary and secondary cal- 
ibrator measurements, which characterizes uncorrected residual 
effects. The total error in the relative flux density measurements 
is typically AS / < S > < 0.5 - 1 %. At /1 9 mm, the increasing 
influence of the atmosphere (and larger pointing and focus er- 
rors) limit our calibration accuracy and the measurement error is 
increased by a factor of up to 2-3. 



2.2. Effelsberg: Total Intensity 

The flux density measurements at the 100 m radio telescope 
of the MPIfR were performed with three multi-horn receivers 
mounted in the secondary focus (see Table |2] for details). The 
4.85 and 10.5 GHz systems are multi-feed heterodyne receivers 
designed as 'software-beam-switching systems'. Both circular 
polarizations (LCP, RCP) are fed into polarimeters providing 
the polarization information (Stokes Q and U) simultaneously. 
In contrast, the 32 GHz receiver is designed as a correlation re- 
ceiver ('hardware-beam-switch') and hence provides only total 
intensity information. 

The target source 07 16+7 14 and all the calibrators are point- 
like within the telescope beam and sufficiently strong in the ob- 
served frequency range. This facilitates flux density measure- 



2.3. Effelsberg: Linear Polarization 

The polarization analysis at 4.85 and 10.45 GHz invokes cor- 
rection of the measured Stokes Q/U channels for instrumental 
polarization, unequal receiver gains, cross-talk and depolariza- 
tion. Since the measurements are done in the horizontal coor- 
dinate system (Alt-AZ), the parallactic angle rotation has to be 
taken into account. For this purpos e we follow the po larization 
analysis using the rnatrix-method of|Turlo et aP (1 19851) (see also 
lOuirrenbach et allll989l:[Kraus et al.Ll2003L for details). We ob- 
tained Stokes Q- and U-amphtudes by fitting Gaussian profiles 
to the pre-averaged sub-scans of each scan for the two orthogo- 
nal and 90 degree phase shifted receiver channels. The 'true' flux 



The results obtained for these sources will be discussed elsewhere. 
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Table 1. The participating radio observatories, their observing wavelengths/frequency, dates and total observing time Tobs- 



Radio telescope & Instittite 


Location 


Aoi,s [mm] 


Vobs LGHzJ 


dates 


'obs LhrsJ 


WSRT (14x25 m), ASTRON 


Westerbork, NL 


210, 180 


1.4, 2.3 


Nov. 10-11 


8 


Elfelsberg (100 m), MPlfR 


Elfelsberg, D 


60, 28, 9 


4.85, 10.45, 32 


Nov. 11-18 


164 


Pico Veleta (30 m), IRAM 


Granada, E 


3.5, 1.3 


86, 229 


Nov. 10-16 


135 


Metsiihovi (14 m), MRO 


Metsahovi, Finland 


8 


37 


Nov. 08-19 


263 


Kitt Peak (12 m), ARO 


Kitt Peak, AZ, USA 


3 


90 


Nov. 14-18 


91 


SMT/HHT(10m), ARO 


Mt. Graham, AZ, USA 


0.87 


345 


Nov. 14-17 


11 


JCMT(15m), JAC 


Mauna Kea, HI, USA 


0.85, 0.45 


345, 666 


Nov. 09-13 


99 



Table 2. Parameters of the EfFelsberg secondary focus receivers. 



System 


4.85 GHz (6 cm) 


10.45 GHz (2.8 cm) 


32 GHz (9 mm) 


Type 


HEMT cooled 


HEMT cooled 


HEMT cooled 


Number of Horns 


2 


4 


3 


Channels 


4 


8 


12 total, 3 X correlation 


System Temp, (zenith) 


27 K 


50 K 


60 K 


Center frequency 


4.85 GHz 


10.45 GHz 


32 GHz 


RF-Filter 


4.6 -5.1 GHz 


10.3- 10.6 GHz 


31 -33 GHz 


IF-Bandwidth 


500 MHz 


300 MHz 


2 GHz 


Polarization 


LHC & RHC 


LHC & RHC 


LHC & RHC 


Calibration 


noise diode 


noise diode 


noise diode 



density vector Smie with the 3 components I, Q and U, which 
describe the source intrinsic polarization, is constructed using a 
3x3-Muller-matrix M, 

S„b,^M-S,n,e^T-P-S,rue ■ (1) 

This allows the parameterization of the parallactic angle rotation 
via the time dependent matrix P and the removal of the afore- 
mentioned systematic effects via inversion of the instrumental 
matrix T (here we assume that Stokes V is negligible small; this 
reduces the matrix dimension to 3 x 3). 

We used the polarization calibrator 3C286 and the highly po- 
larized secondary calibrator 0836-1-710 {p ^ 6% Xo determine 
the elements of the instrumental matrix T via a least- square fit 
procedure to S obs, and by this the instrumental effects. In order 
to improve the accuracy of the determination of the instrumental 
polarization (of the order < 0.5 %), the (non-variable) unpolar- 
ized secondary calibrator 095 1 -1-699 was also used. This proce- 
dure corrects for all (obvious) systematic effects and yields a 
relative measurement accuracy for the polarized flux density of 
^ ~ 5% and Ay = 2-3° for the polarization angle. 

2.4. Observations at higher frequencies 

At the IRAM 30 m telescope, we followed a similar observing 
strategy as in EfFelsberg. In combination with the good telescope 
and receiver performance, a un-precedent precision of the mm- 
flux density light-curve wa s achieved, yi elding an rms-accuracy 
of ~ 1.2 % at 86 GHz (see lAgudo et al.L [2006, for details). At 
229 GHz the calibration accuracy was reduced due to the loss 
of 1 linear feed in the 1.3 mm receiver (mal-function), which 
reduced the measurement accuracy to ~ 16 % rms at 229 GHz. 
This prevents a detailed study of short-term variability at this 
frequency, but still allows the measurement of the overall flux 
density trend. 

The dense in time sampled monitoring at Elfelsberg and Pico 
Veleta was complemented by measurements with some smaller 
telescopes: Kitt Peak (3 mm), SMT/HHT (0.87mm) and JCMT 



Table 3. Total and polarized flux densities and polarization 
angles for the calibrator sources at each observing frequency. 
The quasar 0836-1-714 varies on longer time scales, but on IDV 
timescales is suitable as a secondary calibrator. 



Source 


A 


s, 




X 




[mm] 


[Jy] 


[Jy] 


n 


3C286 


60 


7.49 


0.82 


33.05 




28 


4.45 


0.52 


33.05 




9 


1.84 






3C295 


60 


6.53 


0.01 


71.54 




28 


2.64 


0.04 


-23.61 




9 


0.55 






3C48 


60 


5.50 


0.23 


106.77 




28 


2.59 


0.15 


116.36 




9 


0.80 






NGC7027 


60 


5.47 


0.01 


55.61 




28 


5.94 


0.01 


81.57 




9 


4.57 






0836+710 


60 


2.54 


0.16 


108.13 




28 


2.12 


0.09 


103.37 




9 


1.82 







(0.85 & 0.45 mm). The smaller telescope sizes and the corre- 
spondingly reduced sensitivity and the poorer time sampling 
limit their data quality. For this reason, these data are mainly 
used to extend the spectral coverage towards the sub-mm bands 
and to investigate the basic trend of the variability (see Sect. 3 
and 4). A short summary of the observations with these tele- 
scopes is given below. 

The measurements at Kitt Peak were performed using a 
SIS receiver operating at 90 GHz, while the 345/666 GHz ob- 
servations at the JCMT and the 345 GHz observations at the 
SMT/HHT were performed using bolometer arrays. At the 
SMT/HHT the bolometer consists of 19 individual continuum 
receivers arranged in concentric hexagons, allowing the subtrac- 
tion of atmospheric effects using the off-horns (iKreysa et all 
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10-18 November 2003 A = 6cm 
071 6+7U : Intensity <I> = 1 .594 ± 0.020 Jy 



10-18 November 2003 X = 9mm 
0716+714 : Intensity (I) = 3.528 ± 0.394 Jy 



1.50! 

2954 



1.06 
1 .04 
1.02 
1.00- 
0.98 



2956 2958 2960 2962 

0836+710 : Intensity <I) = 2.542 ± 0.007 Jy 




2958 2960 
J.D. - 2450000 



Fig. 1. Flux density variability of 0716+714 (top panel) at 
60 mm wavelength. For comparison, the secondary calibrator 
0836+710, which shows residual peak-to-peak variations of 
~l% (mo = 0.3 %), is displayed (bottom panel). 



10-18 November 2003 X = 2.8cm 
0716+714 : Intensity <I> = 2.324 ± 0.127 Jy 



2956 

0836+710 : Intensity 



lit 



2958 2960 

<I) = 2.1 19 ± 0.006 Jy 




2958 
J.D. - 2450000 



Fig. 2. Flux density variability of 0716+714 (top panel) at 
28 mm wavelength. The secondary calibrator 0836+710 shows 
residual peak-to-peak variations of ~ 1 % (mo = 0.6 %). 



l2002h . At the JCMT, the continuum array receiver SCUBA was 
used, which has two hexagonal ar rays of bolometric d etectors 
with 37 and 91 horns, respectively (iHolland et al.LlT999h . 

At these three telescopes, the flux density of 0716+714 was 
measured using the on-source-off'-source observing technique 
with the number of on/off subscans matched to the source bright- 
ness at the given frequency. At Kitt Peak and at the SMT/HHT 



^0.9 
^0.8 



2956 

0836+710 ; Intensity 



2958 2960 2962 

(I) = 1.824 ± 0.035 Jy 




Fig. 3. Flux density variability of 0716+714 (top panel) at 9 mm 
wavelength. The secondary calibrator 0836+710 shows residual 
peak-to-peak variations of ~ 8 % (mo - 2.2 %). 



the observing strategy was similar to the one used at Efifelsberg 
and Pico Veleta (see Section 12.11 and 12.21 ). consisting of fre- 
quent observations of primary and secondary calibrators (ultra- 
compact HII regions, planetary nebulae, planets, some AGN). 
Sky dips were done to allow for atmospheric opacity correction 
and the planets were used for the absolute flux density calibra- 
tion (counts-to-Jansky conversion). In all cases, the final mea- 
surement uncertainties are relatively high, up to 15-20%, pre- 
venting us from studying the short-term variability pattern of 
0716+714 at the highest frequencies. We therefore decided to 
improve the signal-to-noise ratio by averaging the individual 
measurements in 24 hr bins. This enables the characterization 
of the main flux density evolution at sub-mm wavelengths over 
the full observing period with a typical measurement accuracy 
of 5-10%. 



3. Data Analysis and Results 

3.1. Total intensity 

In Fig. [U |2] and |3] we plot the flux density measurements of 
0716+714 versus time at 4.85, 10.45 and 32 GHz. In Fig.|4]and 
|5] we report the measurements obtained with the IRAM 30 m- 
telescope at 86 and 229 GHz already published by I Agudo et al.l 
(2006). For a direct comparison, the residual variability of the 
secondary calibrator 0836+714 is shown at each band. 

The light curves given in Fig. [1] to |4] show 0716+714 to 
be strongly variable, when compared to the stationary sec- 
ondary calibrator At /1 60 mm (Fig.[T]|, 0716+714 exhibits low- 
amplitude variations with a slow (~ 2day) total flux density in- 
crease after J.D. = 2452957.6 (peak-to-peak variability ampli- 
tude > 4%). 

The light curves obtained at 28, 9 and 3.5 mm wavelengths 
(Fig. |2]to|4| are dominated by a much stronger mono tonic in- 
crease over a time range of several days between November 
10 - 18. In particular, we note an increasing amplitude of the 
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10-18 November 2003 A = 3mm 
0716+7U : Intensity <I> = 4.499 ± 0.513 Jy 



2954 2956 2958 

0836+710 : Intensity <I) = 



2960 

1.486 ± 0.021 Jy 




2958 2960 
J.D. - 2450000 



Fig. 4. Flux density variability of 0716+714 (top panel) at 
3.5 mm wavelength. The secondary calibrator 0836+710 shows 
residual peak-to-peak variations of ~ 6 % (mo -1.2 %). 



10-18 November 2003 A = 1.3mm 

0716+714 : Intensity <I> = 3.506 ± 0.448 Jy 
6 

5 




2954 2956 2958 2960 

0836+710 : Intensity <I) = 0.830 ± 0.123 Jy 




Fig. 5. Flux density variability of 0716+714 (top panel) at 
1.3 mm wavelength. The secondary calibrator 0836+710 shows 
residual peak-to-peak variations of ~ 65 % (mo - 16.0 %) due to 
a receiver mal-function (see text). 



variability towards higher frequencies (/i28mm: ~16% be- 
tween J.D. =2452956.7 and J.D. =2452961.7; /19mm: -25% 
between J.D. = 2452955.5 and J.D. = 2452958.8). The most dra- 
matic increase is seen at A3 mm (Fig. |4]l. Here, the total flux 
density shows a linear increase of more than 35 % between 
J.D. = 2452954.4 and J.D. = 2452958.7. At A 1.3 mm (Fig. ID, a 
precise characterization of the temporal behavior is not possible 
due to the larger measurement errors. However, a linear fit to 



this data set indicates an increase with a rate of change of 6 % 
per day (see' A^udo et al.L[200 6', for details). 

In order to characterize the multi-frequency variability prop- 
erties more quantitatively, each data set was investigated by 
means of a statistical variability analysis based on the following 
steps: (i) a ;^'^-test for the presence of variability, (ii) the mea- 
surement of the variability strength (modulation index m and 
noise-bias corrected variability amplitude Y), and (iii) the de- 
termination of the characteristic variab ility time scales . Thes e 
met hods are descr i bed in mor e detail by Heeschen et al.l ( Il987h . 
Oui rrenbach et al.1 (1 19921) and lKraus et all (l2003b . In the foUow- 
ing we will use this formalism, and refer to Appendix |A] for 
more details on the definitions of m, and Y. The determi- 



nation of the variabiUty time scales will be presented in Sect. 13.31 
and AppendixlB] 

In Table 4 we summarize the results of the variability anal- 
ysis of the total flux density measurements for 0716+714 (and 
the observed secondary calibrators) at the five observing wave- 
lengths. Here, for each source and each wavelength, the number 
of data points N, the mean flux density < S > and the modula- 
tion index m is given. The (weighted) mean modulation index mo 
of the secondary calibrators is shown in the header of Table 4 for 
each frequency. The corresponding values of mo = 0.3, 0.6, 2.2 
and 1.2% obtained at A = 60, 28, 9, flnt/3mm demonstrate the 
small residual scatter in the secondary calibrator data and thus 
the good overall calibration accuracy. 

The results of the;^'^-test for each data set are given in terms 
of the reduced and the corresponding value at which a 99.9 % 
significance level for variability is reached (xlgg^)- The vari- 
ability amplitude Y is only calculated for those sources which 
according to the ;(f^-test, showed significant variability. It is ob- 
vious from Table 4, that only 0716+714 shows significant vari- 
ability at A 60, A28,A9 and A 3 mm. 

3.2. Polarization 

In Fig.|6]we show the polarization variability of 0716+714 mea- 
sured with the 100 m telescope at /i60mm and /1 28 mm. Here, 
the time evolution of the total flux density / is displayed together 
with the polarized intensity S p and the polarization angle For 
a detailed discus sion of the polariza tion data obtained at Pico 
Veleta we refer to lAgudo et al.l(l2006h . Due to the marginal (2cr) 
detection of polarization variations at 3 mm wavelength, we do 
not consider this in the further analysis. 

A first inspection of the linear polarization data shown in Fig. 
|6] reveals for 0716+714 strong and rapid variability also in the 
polarized intensity and the polarization angle. In order to charac- 
terize the variability behavior in polarization, a variability anal- 
ysis similar to that described for the total intensity in Sect. 13.11 
was performed (see AppendixlAli. The results of the correspond- 
ing quantities {mp, Yp) for the polarization S p are summarized in 
Table 5. A;t'^-test was performed to determine the significance 
of the polarization variability. The corresponding values of 
andxlg 9<jj, for S p andx are also included in the Table 5. 

The statistical variability test shows significant and strong 
variability of 0716+714 also in polarized flux and polarization 
angle. At /160 mm, the polarized intensity displays a 10 - 15 
times higher variability amplitude than the total intensity varia- 
tions. The pronounced decrease of the polarized flux occurring 
after J.D. =2452956.5 (November 13), seems to anti-correlate 
with the increase of the total flux density observed during the 
same time interval. The polarization angle, however, changed 
nearly monotonically by about 20° over the whole observing pe- 
riod. The linear polarization data obtained at A 28 mm exhibit 
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Table 4. Results of the variability analysis for total intensity. N denotes the total number of measurements, < S i > the mean flux 
density with the rms-value os , in the modulation index, Y the variability amplitude, the reduced and;t'gg the corresponding 
value for a 99.9 % significance level of variability (see AppendixlAlfor definitions). 



S5 0716+714 



A 


N 


<5/ > 


as 


m 


Y 


xi 


^99.9% 


[mm] 




[Jy] 


Uy] 


[%] 


[%] 






60 


100 


1.594 


0.020 


1.28 


3.72 


7.12 


1.50 


28 


98 


2.324 


0.127 


5.45 


16.26 


43.87 


1.50 


9 


79 


3.528 


0.391 


11.09 


32.73 


6.38 


1.57 


3.5 


74 


4.502 


0.516 


11.46 


34.19 


55.30 


1.59 


1.3 


68 


3.506 


0.448 


12.77 




0.50 


1.62 



Secondary calibrators 



i60mm m,,o = 0.30%, AlSmm m;,o = 0.55 %, /19mm m,.o = 2.18% 
/13.5mm nijQ = 1.20%, /i 1.3mm m/o = 16.0% 



Source 


A 


N 


<5, > 


as 


m 


Y 


xi 


XggM 




[mm] 




[Jy] 


[Jy] 


[%] 


[%] 






0212+735 


3.5 


33 


1.101 


0.016 


1.45 




0.31 


1.95 




1.3 


15 


0.712 


0.094 


13.14 




0.12 


2.58 


0633+734 


3.5 


25 


0.942 


0.018 


1.90 




0.40 


2.13 




1.3 


15 


0.621 


0.085 


13.72 




0.11 


2.58 


0633+599 


28 


56 


0.677 


0.004 


0.56 




0.46 


1.69 




9 


51 


0.639 


0.015 


2.42 




0.28 


1.73 


0800+618 


60 


94 


1.385 


0.004 


0.31 




0.42 


1.52 




28 


96 


1.130 


0.007 


0.61 




0.54 


1.51 


0835+583 


28 


57 


0.239 


0.002 


0.71 




0.64 


1.69 


0836+710 


60 


98 


2.542 


0.007 


0.28 




0.34 


1.50 




28 


99 


2.119 


0.006 


0.30 




0.13 


1.50 




9 


75 


1.824 


0.035 


1.93 




0.18 


1.59 




3.5 


42 


1.486 


0.020 


1.36 




0.14 


1.82 




1.3 


31 


0.830 


0.123 


14.77 




0.28 


1.99 


0951+699 


60 


48 


3.298 


0.008 


0.25 




0.27 


1.75 




28 


46 


1.812 


0.010 


0.58 




0.51 


1.78 


1803+784 


3.5 


32 


1.334 


0.016 


1.19 




0.09 


1.97 




1.3 


20 


0.940 


0.235 


25.04 




0.78 


2.31 


1928+738 


3.5 


32 


1.866 


0.024 


1.31 




0.18 


1.97 




1.3 


20 


0.904 


0.177 


19.54 




0.44 


2.31 



an overall variability amplitude lower than that of the A 60mm 
data. At A 28 mm, the relative strength of the variations seen in 
the polarized and total flux are comparable, but the polarization 
variations appear faster. We further note that the amplitude of the 
polarization angle variations at A 28 mm are about a factor of 4 
lower than at A 60 mm. Similar as at A 60 mm, the polarized flux 
density variations appear to be anti-correlated with the total in- 
tensity also at A 28 mm. This will be further investigated in Sect. 

EH 

3.3. Variability cliaracteristics 

3.3.1 . Variability time scales for total intensity and 
polarization 

To further quantify the variability behavior of 07 16+7 14, charac- 
teristic variability time scales have been extracted from the light 
curves shown in Fig. [T]-[6] The light curves are dominated by a 
quasi-monotonic flux density increase over several days (and a 
decrease in the polarization), with occasional super-imposed less 
pronounced but more rapid variations. The pronounced, quasi- 
periodic rapid (< 1 day time scale) variability, which is usu- 
ally seen in this a i id oth er IDV sources (e.g. Wagner & Witzel, 
ll995l:lKraus et am2003l) . is not really obvious in this campaign. 
Therefore, a precise determination of a characteristic variability 



time scale (in a statistical sense) is more difficult. We will there- 
fore use a number of different methods to determine the variabil- 
ity time scales: structure function (SF), auto-correlation function 
(ACF) and a minimum-maximum method. These methods and 
the corresponding error calculations are described in more detail 
in AppendixiBland the results will be presented in the following. 
We confine our variability analysis for 0716+714 to the most 
accurately measured data sets and thus exclude the light curve 
obtained at 1.3 mm wavelength. In Table |6] we summarize the 
derived variability time scales and their measurement errors for 
each of the different data sets. 

In Fig.|2]we show examples of calculated structure functions 
for the total and polarized flux density. We note that the error 
of the derived time scales from the SF and ACF analysis shown 
in Table |6] is usually high - values of up to 30% were obtained. 
The uncertainty of the statistical parameters which quantify the 
variations in a time series critically depends on the duration of 
the observation, the sampling interval and the number of signifi- 
cant flux density changes (variability cycles) observed. The light 
curves of 0716+714 (Fig.[T]-|6]l often show only a monotonic in- 
crease and a well defined variability time scale is not observed. 
This allows only to derive lower limits to the variability time 
scales (see AppendixlBli. In the other cases of better defined vari- 
ability time scales, the measured accuracy is a few 10%. Here 
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November 2003 A - 6cm 

0716+714 : Intensity (I) ^ 1-594 ± 0.020 Jy 
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Fig. 6. Flux density variability of 0716+714 obtained at 60 mm (left panel) and 28 mm wavelength (right panel) in total intensity 
(top), polarized intensity (middle), and polarization angle (bottom), respectively. 



Table 5. Results of the variability analysis for polarization. N denotes the total number of measurements, < S p > the mean polarized 
intensity, oSp its rms-value, nip the modulation index, Yp the variability amplitude, < x > the mean polarization angle, cr^ its rms- 
value. The reduced x^- is given for polarized intensity and polarization angle, xlg g% denotes the corresponding 99.9 % significance 
level for variability. See AppendixlAlfor definitions. 



S5 0716+714 



A 

[mm] 


N 


<5p > 

[Jy] 


[Jy] 


mp 
[%] 


Yp 
[%] 


xi 


<X > 

n 


n 


xi 


^99.9% 


60 


98 


0.067 


0.012 


18.14 


54.42 


12.86 


4.40 


4.79 


77.25 


1.50 


28 


98 


0.188 


0.009 


4.78 


14.33 


0.35 


24.92 


2.03 


10.85 


1.50 



Polarization calibrators 



A 60mm 


mp,o = 


1.8% 


(Tj,,o = 0.6° 
















A 28mm 


iripQ = 


2.0% 
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A 
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N 
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[%] 


Yp 
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xi 
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[°] 
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xi 


'^99.9% 


3C286 


60 


16 


0.826 0.009 


1.08 




0.28 


33.02 


0.19 


0.31 


2.51 




28 


15 


0.521 0.008 


1.45 




0.51 


32.95 


0.16 


0.09 


2.58 


0836+710 


60 


95 


0.156 0.002 


1.79 




0.51 


108.17 


0.29 


0.76 


1.51 




28 


97 


0.085 0.002 


1.97 




0.84 


103.56 


0.57 


1.08 


1.51 


0951+699 


60 


48 


0.003 0.002 


















28 


42 


0.002 0.001 

















the values obtained by the three diff'erent analysis methods are 
in good agreement. 

The prominent and monotonic flux density increase seen at 
60, 28, 9 and 3 mm wavelengths appears on similar time scales. 
The SF analysis yields lower limits of 3.7 to 4.2 days (ACF: 3.0- 
3.7 days). Due to the lack of pronounced saturation levels in the 



SFs on time scales < 2 days (see Fig.|7]l, we identify the observed 
inter-day variability of 0716+714 in total intensity as IDV of 
type-I according to iHeeschen et alj (Il987h . A significant com- 
ponent of faster variabiUty - with a time scale of < 1 day - was 
found only at 60mm wavelength. The increasing amplitudes of 
the SFs shown in Fig.|7](left), nicely demonstrate the increasing 
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Structure Functions Structure functions 

lolali,ile.i.s,ly polarised intensity 




time lag [days] 

Fig. 7. Structure functions of 0716+714 plotted versus time-lag. Different symbols are used for different wavelengths. The panel on 
the left shows the SF for total intensity, on the right the SF for polarized intensity is shown. 



Table 6. Variability time scales and lower limits to the variability brightness temperatures for total intensity and polarization in 
0716+714. The variability brightness temperatures were derived using Eq. Q. 
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60 


>6.8 
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28 
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+0.1/-0.2 




3.1 


+0.2/--0.2 











strength of the observed variations towards higher frequencies. 
The value of the SF at small time lags, t — » 0, characterizes the 
noise level present in the time series. We note that for small time 
lags (t < 0.7 days), the value of the A9mm SF in Fig. |7] (top) is 
larger than the A 3 mm SF. This can be interpreted with a reduced 
short time stability of the 9 mm signal, which is affected by the 
stronger atmospheric influence at the EfFelsberg site compared 
to the high-altitude Pico Veleta site. 



3.3.2. Variability time scales and brightness temperatures 

Assuming that the observed rapid variability would be source 
intrinsic (see Sect. 14.21 ), the variability time scales summarized 
in Table|6]imply very compact emission regions, and hen ce very 
high intrinsic brightness temperatures. Following Marsche r et al.l 
(Il979l) , and taking an isotropically expanding light-sphere with- 
out preferred direction into account, the light travel time argu- 
ment implies a diameter d < 2 ■ cA f of the emitting region after 
the time interval At [s] of expansion. We then obtain for the di- 
ameter of the variable emission region [mas]: 



In linear polarization, the source varies at A 60 mm on time 
scales comparable to those seen in total intensity, whereas at 
A 28 mm, two variablity time scales are seen: one comparable 
to the time scale of the total intensity variations (~ 4 days) and 
a significantly faster component, showing f^^ - 1.3 days. The 
polarization angle changes at A 60 mm display a continuous de- 
crease over the seven days and our analysis yields a lower limit 
of tvar > 6.8 days. In contrast, the variations of the polarization 
angle at A 28 mm appear faster by a factor of about two. 



61 = 0.13— (5(1 +z) , (2) 

where z is the redshift of the source, 6 the Doppler factor, ?,,£»• 
the logarithmic variability time scale in years (see Appendix iBl 
eq. B5) and di the luminosity distance in Gpc. Consequently, we 
find for the brightness temperature Tb [K] of a sta tionary compo- 
nent with a Gaussian brightness distribution (e.g. lKovalev et all 
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Fig. 8. Examples of cross-correlation functions between total intensity light curves obtained at the different frequencies. Here, we 
show combinations with respect to the data at 60 mm (see top right labels in each plot). The CCPD for the peak of the DCF is shown 
in the right-hand panels (for details see text). In the bottom right panel, we plot the time lags versus frequency, relative to the light 
curve at 60 mm. Note the trend of increasing (negative) time lag with increasing frequency. 



I2005h and flux density Sa [Jy]: 



Tb = 8.47 X lO'*- 5, 



A dr. 



(3) 



According to Eq. (O we calculated the brightness temperatures 
for the 'most reliably determined' variability time scales derived 
in the previous section (see also Appendix |B]). Here and in the 
following we adopt a cosmological distance corresponding to 
z = 0.3, which is a lower limit for the redshift of the source 
(see Sect.[T), and yields a luminosity distance of di = 1-51 Gpc 
assuming a flat universe w ith i/n - 72 k r n s~ ' M pc~ ' , Q.m = 0.3, 
Qa = 0.7 and = dSpergel et al.L l2003l) . In Table |6l we 
summarize the calculated values of Tg for each observing band 
and for total and polarized intensity, respectively. 

For the slow variability of type 1 seen in the cm-bands, the 
brightness temperatures range between 4.010''* K and 4.9-10'^ K, 
whereas in the mm-bands, we obtain lower values ranging be- 
tween 6.5 ■ 10'^^ K and 4.2 ■ 10"* K. The faster type 11 variability 
component seen at A 60 mm, leads to a brightness temperature of 
up to 8. 1 ■ 10"' K, which is at least one order of magnitude higher 
than the previous estimates. The variations observed in polar- 
ized intensity reveal brightness temperatures ranging between 
3.2 ■ 10''*Kand2.1 ■ 10'^ K, similar to those obtained for to- 



tal intensity in the cm-band. We note that in our calculations we 
used conservative estimates of the variability time scales, which 
yield hard lower limits to Tb- Depending on the assumptions 
on the geometry and isotropy of the emitting region, the actual 
brightness temperatures may be higher by a factor of up to 6 
(e.g. assuming for the emission region a uniform flat disk of size 
d ^ cAt as in the case of a shock). We also note that a larger 
cosmological distance of 0716+714 would further increase Tb- 



3-4- Cross correlations 

To further quantify the close similarities seen in the light curves 
of 0716+714 across the observing bands, we computed the dis- 
crete cross-correlation function (DCF) between different fre- 
quencies allowing to search for correlations and possible time 
lags. W e f ollowed the method describe d by Edelson & KroliS 
(1 19881) and iHufnagel & BregmanI (Il992h for unevenly sampled 
data, and calculated the DCF and the position of its maximum 
by using the centroid Tc of the DCF, given by 



ZiDCFi 



(4) 
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In order to obtain statistically meaningful values for the cross- 
correlation time lags and their related uncertainties, we per- 
formed Monte Carlo simulations, fo l lowin g the methods de- 
scribe d in detail by iPeterson et alj (Il998l) and iRaiteri et"aLl 
(I2003h . Here, we take into account the influence of both uneven 
sampling and flux density errors. This was done using random 
subsets of the two data sets. In addition, random Gaussian fluc- 
tuations constrained by the measurement errors were added to 
the flux densities. In each simulation we determined the centroid 
Tc of the DCF peak. After running 1000 simulations, we obtained 
a cross-correlation peak distribution (CCPD), which is shown in 
the right panel of Fig. [8] This technique yields a reliable measure 
of the uncertainties in the estimated time lags whereas the val- 
ues, cornputed_directl};jrom the CCPDs, correspond to 1 cr errors 
(see lPeterson et all [199 8!). 

This procedure was applied to each possible frequency 
combination of the total intensity data (60/28 mm, 60/9 mm, 
60/3 mm, 28/9 mm, 28/3 mm, 9/3 mm). Fig. [8] shows some ex- 
amples of the DCFs and CCPDs relative to the /1 60 mm data. 
Our analysis confirms the existence of a significant correlation 
across all observed radio-bands. This enables the determination 
of time lags, which are T6o/28 = ~0-3!^o 3 days between the 60 mm 
and the 28 mm data, T(,o/9 = -2.1 '^^^ days between the 60 mm and 
the 9 mm data, and r6o/3 = -2.7^g 4 days between the 60 mm and 
the 3 mm data. This systematic trend is also seen in the time lags 
relative to the 28 mm data, with T2S/9 = -1.9;^g^days and T28/3 - 
-2.4;^y 2 days, respectively. In Fig. [8] (bottom right) we plot the 
Tfio/i time lags for the 3 observing bands (/ = 28, 9, 3 mm) versus 
observing frequency. A systematic trend of increasing (negative) 
time lag towards higher frequencies is evident. This shows that 
the observed variations (the observed flux density increase) oc- 
cur first at higher frequencies and then propagate through the ra- 
dio spectrum towards lower frequencies. The time delay between 
the two most separated bands (A 3 mm and A 60 mm) is about 
2.5 days. This time-lag behavior together with the observed in- 
creasing variability amplitudes will be discussed in Sect. 4. 

In order to investigate the possible anti-correlation of the to- 
tal intensity and polarization variations mentioned earlier, we 
performed a cross-correlation analysis (DCF) also between to- 
tal and polarized intensity, and between polarized intensity and 
polarization angle for the 60 and 28 mm wavelengths data. The 
results are shown in Fig. |9] For both radio bands we confirm an 
anti-correlation between the total and polarized intensity with 
a trend of the total intensity leading the polarized intensity by 
0.5 days at 60 mm but no obvious time lag at 28 mm. Formally, 
we also find an anti-correlation between the polarized intensity 
and polarization angle at 28 mm wavelength. 

3.5. Broad-band spectra 

In Fig. [TOl we show the combined variability data obtained for 
0716-1-714. This multi-frequency data set allows to study the 
spectral evolution and variability on a daily basis. We construct 
(quasi-) simultaneous radio spectra using daily averages from 
the data of all participating radio observatories: Effelsberg, Pico 
Veleta, SMT/HHT, JCMT, KP and WSRT (see Table 1). The 
spectral coverage ranges from 1 .4 GHz to 666 GHz. At 3 mm 
wavelength we combined the data from IRAM and KP, to ex- 
tend the time coverage up to November 17. At 0.8 and 0.85 mm, 
we used the data from JCMT and the SMT/HHT, which unfor- 
tunately do not overlap in time. The time evolution of the cm- to 
sub-mm spectrum of 0716+714 over the seven observing days 
(November 11-17; hereafter referred to as days 1 to 7) is pre- 



sented in Fig. [14] Due to the diff'erent duration of the observa- 
tions at the difi'erent telescopes, the maximum frequency cover- 
age was obtained only on day 1 (1 .4-666 GHz), whereas for days 
6 and 7 the frequency coverage is reduced to 4.85-345 GHz. 

0716+714 shows a very inverted radio spectrum over the 
whole observing period, with small, but significant brightness 
variations occurring near the spectral turnover during the first 
five observing days. A pronounced spectral maximum is seen 
near 90 GHz. Linear fits to the daily spectra yield an averaged 
spectral slope a,i,icii {Sv~ v") of + 0.39 + 0.02 for the optically 
thick part of the spectrum between 1 .4 and 86 GHz. lAgudo et al.l 
(2006) give for the spectral slope between 86 GHz and 229 GHz 
a typical (ise/aiQ - -0.23+0.10, a characteristic for the transition 
towards optical thin synchrotron emission. With the broader fre- 
quency coverage shown in Fig. [14] we obtain an average value of 
(ithin - - 0.30 ±0.05. We note that a more detailed analysis of the 
high frequency data from JCMT indicates that the calibration of 
the 666 GHz data is uncertain and the measured fluxes come out 
very low. For this reason we did not include the 666 GHz data in 
the spectral fitting. 

The spectral slopes obtained from the daily fits suggest slight 
changes with a moderate spectral hardening in the optically thick 
part (from a/hick - +0.37 to a, hick - +0.41 between days 3 and 
4) and a spectral steepening in the optical thin part (from a,/„„ = 
-0.25 to a,hi„ = -0.38 between days 1 and 3). For the turn- 
over frequency v,„ we made parabolic fits to the daily spectra, 
which suggest a small shift of v„, of about 10 GHz towards lower 
frequencies. However, these changes in a and v„, appear not to 
be statistically significant within our measurement accuracy and 
the change in frequency coverage of our spectra. A more detailed 
discussion of the spectrum and its variability will be presented 
in Sect. [43] 

4. Discussion 

4. 1 . Variability ciiaracteristics 

4.1 .1 . Variability in tine cm- to mm-regime 

During this campaign 0716+714 exhibits very different variabil- 
ity characteristics when compared to previous IDV observations 
in the cm-regime (e.g. lOuirrenbach et al.L Il992t iKraus et all 
120031) . Usuafly the source shows rapid IDV on time scales of 
< 0.5-1.5 days (in intensity and polarization). In this observa- 
tion, however, it showed slow inter-day variability on a > 3-4 day 
time scale (total intensity). When 0716+714 exhibits rapid IDV, 
the variability amplitud es often decrea ses between 5 and 10 GHz 
(see Fig. 1 of Krichbau m et al.L |2002|) . A variability amplitude 
constant or increasing with frequency was observed only oc- 
casionally. It is therefore remarkable that this new observations 
show such a systematic and strong increase of the variability am- 
plitudes (from F = 3.7 % at 5 GHz to 7 = 34 % at 90 GHz). 

A direct comparison with /1 60 mm Effelsberg IDV light 
curves obtained earlier (April 2002) and 8 month later (July 
2004) reveal the 'classical' IDV behavior, with variability time 
scales of 0.5-1 days and modulation indices of 2. 1-3.4 %. These 
data sets are shown in Fig. [TT] for comparison. This suggests 
possible changes of the variability characteristics over time 
scales of weeks to months. We note that a change of the vari- 
ability time scale fror n ~ 1 to ~ 7 days was also observed 
in 0716+714 in 1990 (lOuirrenbach et al.1 Il99ll; IWagner et all 
1996). Other prominent examples of changes of the variability 
mode or episodic ID V are PKS 04 05-385 (Kedziora-Chudczeri 
2006) and 0917+624 (iKraus et al.L[T999HFuhrmann et alifeOOa 
Fuhrmann et al., in prep.). 
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Fig. 10. Summary of the total flux density variations in 0716+714 obtained at the various cm- to sub-mm bands during the time of 
the core campaign. The contribution from each observatory is indicated in the top right insert and in Table [1] (see text for details). 
For better illustration, the light curve at /1 1.3 mm (see Fig.|5]), which is affected by larger measurement errors, is not shown. 



The good frequency coverage of the observations presented 
here provides for the first time a possibility to study with dense 
time sampling the short-term variations from the cm- up to the 
mm-regime. Here, the correlated variability across all bands, 
the observed frequency dependence of the variability ampli- 
tudes, and the observed time lag with variations at higher fre- 
quencies appearing earlier (Sect. 13.41 Fig. [8]i, argue in favor 
of a source-intrinsic origin of the observed variations. Such 
'canonical' variability behavior is usually observed in AGN 
and other compact radio sources over longer time intervals 
(weeks to years) and is commonly explained by synchrotron- 
cooling and adiabatic expansion of a flaring co mponent or 
a shock (e.g. van der Laan, 1966; K ellermann & Pa ulinv-Totm 
ligeS ^.'Marsc her & Gear,, J985). In this framework. ' Raiteri et all 
(1200 3) found similar variability characteristics for the long-term, 
multi-frequency radio light curves of 0716+714. A possible con- 
tribution of scintillation effects and thus extrinsic origin will be 
discussed in Sect. 14. 21 A co mparison with the simu ltaneous opti- 
cal R-band data presented bv lOstorero et alj ( l2006l) will be given 
in Sect. |4T2l 

The variability brightness temperatures derived for the total 
intensity data sets (Table |6] Sect. 13.3.21 ) reveal lower limits of 
4 - 8 ■ 10"' K for the faster variability component observed at 
60 mm wavelength, and 10'^ ** - 10'^ K for the correlated cm- 
/mm-wavelength flux density increase on inter-day time scales. 
We notice a systematic decrease of Tb towards higher frequen- 
cies. Since Tb was derived using the light travel time argument, 
which implies upper limits to the size of the region of variable 
emission, it remains uncertain if the observed frequency depen- 
dence of the lower limits of Tb also reflects a similar frequency 
dependence of the source-intrinsic brightness temperature. The 
apparent bri ghtness temperatures, however, exc eed the IC limit 
of ~ 10'2 K (iKeHermann & Paulinv-TotR Il969t) by at least ~ 2 
orders of magnitude, even at mm- wavelengths. If the excessive 
Tb is the result of relativistic aberration, the source radiation 
must be strongly Doppler-boosted in all of the observed wave- 



bands. The resulting Doppler-factors will be discussed in Sect. 
El 

From the analysis presented in Sect. 3, we deduce a com- 
plex behavior of the polarization variability. We find (i) more 
pronounced variability amplitudes in polarization than in total 
intensity (a factor of ~ 15 at /1 60 mm), (ii) more pronounced po- 
larization variability (in P and x) at ^ 60 mm than at A 28 mm, 
(iii) significantly faster variability at A 28 mm, and (iv) a clear 
anti-correlation between total and polarized flux density at both 
observing bands. Such a complex variability beh avior in po- 
larization is often observed in IDV sources (e.g. iKraus et all 
|2003), and is interpreted by a multi-component sub-structure of 
the emitting region(s). The superposition of individually vary- 
ing and polarized sub-components (characterized by their mis- 
aligned polarization vectors) could in principle reproduce the ob- 
served polarization variations. At present it is unclear, whether 
the necessary variability of the individual sub-components is 
source-intrinsic (e.g. J^ian, 1993)) , sole l y cau sed by interstel- 
lar scintillation (e.g. lRickett et al., 1995, 2002), or results from 
a mixture of these efifects 40ian et al .. 2002)- The recent de- 
tection of polarization IDV in the VLBI core of 0716+714 
clearly places the origin of IDV in the unresolved VLBI core- 
region, and strongly supports the idea of multiple emb edded 
sub-components of < lOOjuas size (see B ach et al.L 12006). It is 
unclear, however, whether the different time lags seen between 
the total and polarized intensity at /1 60 mm and /1 28 mm have a 
physical meaning. Source-intrinsic opacity effects and a possible 
spatial displacement of the polarized sub-components relative to 
the total intensity component, as often obser ved for ind i vidual 
jet components with polarization VLBI (e.g. i Bach etalil2006l 
and references therein), might provide a simple explanation. 

4.1.2. Long-term variability and broad band cliaracteristics 

During the variability campaign in November 2003, 0716+714 
was observed in a particular active phase, with a high radio- 
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Fig. 9. Cross-correlation functions between total intensity and 
polarization at 60 mm and 28 mm wavelengths. The top panels 
show the DCFs between total and polarized flux density. The 
bottom panel shows the DCF for polarized intensity and polar- 
ization angle at A 28 mm. 



to-optical state. In fact, between September and October 2003 
the source underwent a dramatic and unprecedented outburst in 
the cm- to mm-bands. In Fig. [12] we show the long-term light 
curves obtained between 2001 and 2004 at frequencies ranging 
between 5 and 90 GHz. In addition to the data from this pa- 
per, the figure also includes data from the Michigan monitor- 
ing (A ller & AUer, priv ate comm.), the Metsahovi blazar moni- 
toring (Terasranta et al., 2005) and from IRAM (H. Ungerechts, 
private comm.). Albeit being in the declining phase of the previ- 
ous large outburst (peak in 2003.8), 0716-i-714's flux was rising 
a second time and was particularly bright in November 2003, 
reaching more than 5 Jy at 3 mm wavelength. The overall in- 
verted radio-to-mm spectrum during this time indicates a rela- 
tively high opacity of the source (synchrotron self-absorption), 
which is a possible cause for the absence of the usually more 
pronounced IDV in 0716+714 at the longer cm- wave lengths. 

On VLBI scales 0716+714 frequently shows component 
ejections on time scales of ~ 1-2 years w hich often are pre - 
ceeded by larger flux-density outbursts (e.g. iBach et al.L l2005l) . 
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Fig. 11. IDV fight curves of 0716+714 observed with Efi'elsberg 
in April 2002 (top) and July 2004 (bottom) at A 60 mm. Note the 
different and more typical ('classical') variability pattern, with 
faster (~0.5-l day) and more pronouncedm - 2.1-3.4%) vari- 
ations than those observed in November 2003 (Fig. [1]). In the 
cases shown here, the source was in a fainter radio state than 
in November 2003, with a mean total flux density of 0.69 and 
1.07 Jy in April 2002 and July 2004, respectively. 



Our daily VLBA observations performed during November 1 1- 
16, 2003 (six epochs), however, indicate no strong structural 
changes at 22 and 43 GHz and on the mas-scale. The flux density 
of the VLBI core (a ~ 35 % increase from day 1 to 6 at 43 GHz) 
seems to follow the trend which is seen in the total intensity light 
curves. A more detailed study of the VLBI results will be pre- 
sented in a future paper (Agudo et al., in prep.). 

During the campaign, 0716+714 was recorded in a moder- 
ate level of optical emission in agreement with previous find- 
ings, where no obvious strong correlation between prominent 
radio and optical outbursts have been observed in this source 
(Raiteri et al., 2003). In the optical R-band, 0716+714 showed 
intra-night to inter-day variations. The root mean square vari- 
ability amplit ude during the time of the core campaign is ~ 23 % 
(see Fig. 1 in lOstorero et al.L[2006h . No obvious correlation be- 
tween the optical and radio variability is seen, indicating that 
either the emission regions generating radio and optical vari- 
ability are spatially separated, or the variability in both bands 
has a different physical origin. We note that also no obvious 
simultaneous (or delayed) transition of the variability mode is 
seen bet ween radio and optical, similar to the one reported by 
Ouirrenbach et alj (11991.) . It is remarkable that during this latter 
campaign, the source was in a much fainter radio-to-optical state 
and the radio spectrum (between 5 and 8 GHz) was much less in- 
verted (Wagner et al., 1990; Oian et an,ll996h . The presence (or 
absence) of radio-optical correlations may depend on the source 
opacity, which is directly related to the observed shape of the 
radio spectrum and the state of source activity. 

Furthermore, we note that in November 2003, 0716+714 
was not detected by the INTEGRAL satellite at any of its 
high energy bands (3-200 keV) and only upper limits to the 
hard X-ray/y-ray emission could be obtained (see Section |4.4.5l 
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lower than in mos t earlier experiments (e.g . lOuirrenbach et all 
[1991 iKraus et al.L llool [Fuhrmanni l2004t iBach et al.L 120051) . 
Typically, the variability index peaks around 5 GHz and de- 
creases towards both higher and lower frequencies. This is ex- 
pected for the transition between weak and strong scattering. The 
strength of the variability observed in November 2003, however, 
strongly increases with observing frequency (see Sect. [3] and H] 
and Fig.[T3Tl. (ii) 0716+714 usually varies fast, on time scales a 
factor of ~ 2-10 faster than observed here. The prolongation of 
the variability time scale in November 2003 might be explained 
as a seasonal effect, caused by the annually changing align- 
ment of the velocity vectors betwe en Earth and s cattering screen 
(as ob served for e.g. J1819+3845: lDennett- Thorpe & de BruvlH 
I2002h . Screen velocities of about 25 km s ^ w.r.t. the local stan- 
dard of rest would be sufficient to explain a prolongation in 
November However, the multi-frequency time lag (correlation) 
of the variability pattern, the increase of the modulation index 
from cm- to mm-bands and the low variability amplitude at 
A 6 cm are not in agreement with such a scenario. As purely ge- 
ometric effect, the orbital motion of the Earth's can affect the 
variability time scale, but not the strength of the observed varia- 
tions. We also note that 0716+714 is observed every 4-8 weeks 
in a regular IDV monitoring project performed with the Urumqi 
telescope. So far, no strong e vidence for annual m odulation of 
the variability pattern is seen (iMarchih et al.Ll2008l) . 

At frequencies of a few GHz, the angular size of an ex- 
tragalactic self-absorbed, incoherent synchrotron source is usu- 
ally expected to be larger than the Fresnel scale Op for a scatter- 
ing screen beyond the Solar system (e.g. iBeckert et all |2002|) . 
Consequently, we will use a slab-model for weak ISS assuming 
quenched scattering as described by Beckert et al. (2002). In the 
regime of weak, quenched scattering (6 > Of), the authors give 
the following analytical solution for the modulation index m: 
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Fig. 12. Long-term variability of 0716+714 at 5, 8, 15, 22, 37 
and 90 GHz, showing flux densities observed between 2001.5 
and 2004 with the Effelsberg, Metsahovi, Michigan and Pico 
Veleta radio telescopes. The data sets were taken from Terasranta 
et al. 2005, Aller & Aller (private comm.), H. Ungerechts (pri- 
vate comm.), and this paper Note the unusual source brightness 
during the outburst around 2003.75. The time of the core cam- 
paign (November 8-19, 2003) is marked by dashed vertical lines. 
For better illustration, an enlargement of the variability curves 
for the time range 2003.45 - 2003.95 is shown in the bottom 
panel. 



and lOstorero et al.L l2006h . However 0716+714 was seen by 
INTEGRAL a f ew months later (at 30-60 keV) (April 2004, 
|PianetaUl2005l) . and some flux density variability in the X-ray 
band was noted. 



4.2. Variability versus frequency: comparison with weak 
interstellar scintillation 

Local scattering in the ISM was clearly demonstrated to be the 
origin of the most rapid (sub-hour) variabili ty see n in a few IDV 
sources (e.g. Dennett-Thorpe & de Bruyri, 1200 lb . Further, ISS 
on longer time scales (da ys to weeks) is prese nt in a large num- 
ber of radio sources (e.g. i Rickett et al.L l2006h . Here, we inves- 
tigate a possible contribution of ISS to the observed cm- up to 
mm-variability of 0716+714. 

In November 2003, 0716+714 varied in a different manner 
than in previous IDV observations: (i) The observed variability 
index m at /1 60 mm (1.3 %) is (by a factor ~ 2-4) significantly 



m^(A) 



2(^) A\D6f-^SM-fm 



(5) 



with the distance to the scattering screen D, a circular source 
model with size 6, the scattering measure SM, the power law 
index /3 = 11/3 for a Kolmogorov spectrum of density irregu- 
larities and the function /i (J3) of order unity. Depending on the 
frequency dependence of the source size and for high galactic 
latitude sources (as 0716+714), the amplitude of the ISS induced 
variations should thus strongly decrease towards higher frequen- 
cies. 

Assuming a slab of 1 pc thickness for a thin screen with 
a strength of turbulence C^, similar to the values derived 
from pu lsar scintiUation for enhanced scattering in the Local 
Bubble ( Bhat et al.L 1199 8). we qualitatively compare our multi- 
frequency results of m directly with ISS according to Eq. 
(|5]l. Combined Space-VLBI and single-dish observations of 
0716+714 place the physica l origin of IDV inside the VLBI- 
core region lBach et al.l (l2005h . which contains about 70 % of the 
source's total flux density. Thus, we assume a scintillating com- 
ponent containing 70 % of the observed total flux density, and 
re-scale the observed modulation indices given in Table 13.11 ac- 
cording to »J70% = m/Q.l. In this way, we obtain mj()% = 1.8 % 
at A 60 mm, which can be reasonably reproduced by assuming 
a scintillating source size of 0.25-0.3 mas and and a screen dis- 
tance of 100-200 pc. Taking these values and assuming a linear 
dependence of with frequency, we can compare the observed 
modulation indices at higher frequencies with those expected 
for weak ISS. The results are shown in Fig.[T3] A comparison 
between observations and model strongly suggests that the ISS 
model can not reproduce the increasing variability amplitudes 
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Fig. 13. Variability amplitudes of 0716+714 across the observ- 
ing bands: observed modulation indices versus frequency with 
ISS model predictions for different parameter ranges (screen dis- 
tance and source size at A 60 mm) superimposed. Note the strong 
discrepancy between model calculations and observed values 
(see text). The statistical errors of m were obtained taking into 
account the limited total observing time and sampling interval 
(see Gabanyi et al. 2007). 

observed towards higher frequencies. We obtain similar results 
by systematically changing the model parameters such as screen 
distance, source size and scattering measure. Consequently, we 
conclude that the observed simultaneous and correlated inter-day 
variations of 0716+714 on time scales of several days in the cm- 
to mm-band are not strongly affected by scintillation effects, and 
therefore should be considered as intrinsic to the source. 

If ISS dominates the variability at A 60 mm, however, the low 
variability amplitude and the long variability time scales (com- 
pared to previous observations) must indicate strong changes, 
either in the scattering medium, or in the intrinsic source struc- 
ture on scales comparable to the scattering size of the medium. 
Changes in the scattering properties of the screen can be re- 
lated to either changes in the distance or the strength of turbu- 
lence. Both appear unlikely to happen on timescales as short as 
months. Since 0716+714 was in a flaring state during this ob- 
servations (see Sect. 14. IK changes in the source structure there- 
fore appear more likely. Flux density outbursts are often ac- 
companied by the ejection of new jet components which, for 
instance, would temporarily increase the size of the scintillat- 
ing core region. In such a scenario, the size of the scintillating 
component should have increased by a factor 2-3 (from about 
0.08-0.09 mas to 0.25-0.3 mas) in order to quench m, down to 
the observed level of 1.8 %. The observations performed eight 
months later (July 2004, Fig. fTTT i reveal a higher modulation in- 
dex of mjo% - 3.0%, indicating again a smaller scintillating 
source size of about 0.1 7-0. 19 mas (assuming a screen distance 
of 100-200 pc). Future VLBI studies with highest possible angu- 
lar resolution will be needed to relate component ejection, mo- 
tion and core size evolution with the changes of the IDV pattern 
in 0716+714. 

4.3. Synchrotron spectra and spectral evolution 

In Fig. [14] we show the evolution of the cm- to sub-mm radio 
spectra on a daily basis. All spectra are strongly inverted, peak- 
ing at 4 - 5 Jy near 3 mm wavelength (90 GHz). The results of 
the spectral fits described in Sect. 13.51 demonstrate that during 
the whole core campaign the radio spectrum was optically thick 
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Fig. 14. Simultaneous cm- to sub-mm spectra of 0716+714 de- 
rived from daily averages of the flux-density measurements 
shown in Fig.[TO] The top-panel shows the spectra for days 1-4 
(November 1 1-14, J.D. = 2452955-58), whereas the spectra of 
days 5-7 (November 15-17, J.D. = 2452959-61) are shown be- 
low. The flux density errors are small and when invisible, com- 
parable or smaller than the size of the symbols. 



up to the mm-band with aiinck - + 0.39 + 0.02 and optically thin 
{(iihin - - 0.30 + 0.05) beyond a synchrotron turnover frequency 
Vm somewhere close to 9 GHz. This confirms previous findings 
by Ostor ero et al.l (12006 ). who showed the broad band spectral 
energy distribution (from radio to y-rays) using the data of this 
campaign. The radio spectra shown here (Fig. [I4] | largely dif- 
fer from earlier multi-frequency studies, where the source typ- 
icall y showed a flat and only mildly inver t ed ( 1 Jy) spectrum 



I showed a ilat and only mildly mver t ed (iJy) spectrum 
Chini et al.L [1988; Wa gner et all '1990'; O uirrenbach et all 



Il991h . In contrast, the new spectra are indicative of a high peak- 
ing and self-absorbed flare spectrum. 

Significant changes in the cm- to sub-mm spectra occur 
mainly until day 5, i.e. during the period of the monotonic in- 
crease in total flux density seen across all bands. Here, the spec- 
tral peak flux increases, reaching its maximum on day 5. The 
daily averaged peak flux density (at v^) rises from Ss6CHz - 
4.0 Jy on day 1 to Ss(,chz - 5.1 Jy on day 5. After this date, no 
further changes are evident (see also Fig.fTOb. 

The absence of pronounced variations in the spectral slopes 
and the turnover frequency (achromatic variations) would indi- 
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cate a geometrical origin, i.e. caused by changes in the beam- 
ing factor (e.g. in a helical/precessing iet:'Begelman et aLlll980t 
[Villata & Raiteri, 1999). The possible slight changes in the daily 
spectral slopes (spectral steepening for v < v„„ the spectral soft- 
ening for V > v,„) and marginal shift of the turn-over frequency 
towards lower energies (Av ~ 10 GHz), suggest that the spec- 
tral evolution is not fully achromatic. The limited time cover- 
age (only 7 days) prevents a more detailed analysis of the spec- 
tral evolution. However, the presence of source intrinsic spec- 
tral variations is also supported by the canonical behavior of 
the time lag (discussed in Sect. 13.4b and the related frequency 
dependence of the variability amplitudes, which increase to- 
wards the mm-bands. Therefore, an interpretation within stan- 
dard flare/adiabatic expansion and/or shock- in-jet models ap- 
pears applicable (e.g. ' Marscher&Ge^ll985l) . 

The optically thick spectral index, however, strongly dif- 
fers from the canonical value for a self-absorbed, homoge- 
neous synchrotron component (Sv ~ v^^^). This indicates that 
the flaring component, which dominates the observed radio- 
spectrum, is inhomogeneous and may consists of more homo- 
geneous but smaller sub-components, pea king at slightly differ - 
ent self-absorption frequencies (see also Ostorero et al.L l2006l) . 
This view is further supported by (non-simultaneous) mm-VLBI 
observations of 0716H-714, which indeed show multiple compo- 
nents on the 0.05 -0.1 mas scale (T. Krichbaum, unpublished 
data). 

4.4. Doppler factors and magnetic field 

In this section we will compare the Doppler factors dyar deduced 
from the variability brightness temperatures of Sect. 13.3.21 with 
Doppler factors calculated by other methods. In particular, we 
will compare d^ar with (i) estimates derived from VLBl studies 
of 0716+714 (dvLBi), (ii) the equipartition Doppler factor deq, 
using calculations of the synchrotron and equipartition magnetic 
field, (iii) dmr.eg using the equipartition brightness temperature 
Tgg, and (iv) the inverse Compton Doppler factor die as calcu- 
lated from the simultaneous INTEGRAL observations. 



4.4.1. Doppler factor 5„ar,/c 

The limits of the apparent brightness temperature {7"^'^) sum- 
marized in Table |6] exceed by several orders of magnitude th e 
IC-limit of ~ 10'^ K dKellermann & Pauhny-Tothl Il969h . 
In order to explain the excessive temperatures solely by rela- 
tivistic boosting of the radiation, source intrinsic radiation near 
the IC-limit can be assumed (Jb - Tg'"^), which leads to 
a lower limit to the Doppler factor of the emitting region of 

Svar,ic - (1 + z) '*^r^'^''/10'^. This leads to Doppler factors 

5var,cm > 10-22 in the cm-regime and 5„ir.mm > 5-10 in the mm- 
regime (assuming z = 0.3), which lower the observed bright- 
ness temperatures of Table |6] down to the IC-limit. In the case 
of the faster variability component seen at /160 mm, the situa- 
tion becomes more extreme. Here, Doppler factors dyar.ic >43 
would be needed if, again, sole intrinsic in origin. Although 
high, these limits for dyar.cm and 6yar,mm agree well with recent 
estimates coming from kinematic VLBI studies of the source. 
From 26 epochs o f VLB I data observed between 1992 and 2001, 
iBach et al] (l2005h found atypically high superluminal motion 
and conclude that the minimum Doppler factor is 6vlbi ^ 20-30, 
assuming a redshift of z > 0.3. This confirms the previously re- 
ported high speeds (16-21 c) in 0716+714 ( Jorstad et al., 2001). 
Moreover, using a larger sample of sources. Jorstad et al.i(i2005i) 



observed 6vlbi > 16-21 in 9 of the 13 studied blazar indicating 
that strong Doppler boosting is not a rare phenomenon in com- 
pact VLBI jets. 

We further note that the observed decrease of Tg'''' towards 
higher frequencies suggests that the Doppler factors at cm- and 
mm-wavelengths could be different. This appears not unlikely 
in a model of a stratified and optically thick self-absorbed bent 
jet, where at each wavelength spatially different jet regions are 
probed. Each of these regions may exhibit characteristic Doppler 
factors depending on jet speed and viewing angle. As one looks 
at higher frequencies into regions which are deeper embedded in 
the source, a lower brightness temperature and resulting lower 
6yat- towards the innermost jet regions would then indicate ei- 
ther jet-acceleration (changes of the Lorenz factor along the jet) 
and/or jet bending with changes of the jet orientation away from 
the observer's line of sight for outward oriented motion. 



4.4.2. Magnetic field from synchrotron self-absorption 

Using standard synchrotron expressions, it is possible to con- 
strain the magnetic field B in 0716+714 by means of the spec- 
tral parameters deduced in the previous section. Assuming a ho- 
mogenous synchrotron se lf-absorb ed region, a closed expression 
for B is given by (see e.g. lMarsch er. 1987) 



BsA [G] = 10-' b(a)ef vis 



(6) 



with the tabulated quantity h(a) depending on the optically thin 
spectral index Othin (see Table 1 in lMarschen,ll987 *). the flux den- 
sity S ,„, the source angular size 6 at the synchrotron turnover fre- 
quency v,„ and the Doppler factor S. At v,„ - 86 GHz, the size 
of the emitting region responsible for the observed variations 
can be constrained using recent mm-VLBI measurements of the 
core region of 0716+714 (B ach et al.L 120061: lAgudo et al.Ll2006l) 
which yield 0%(,chz < 0.04 mas. Using the observed spectral pa- 
rameters from above {S%^ghz =4.0 - 5.13 Jy, athin - - 0.30,) we 
calculate a lower limit of the magnetic field Bsa in the range of 
0.07 - 0.1 1 5G for v,„ = 86 GHz without correcting for Doppler 
boosting. 

4.4.3. Magnetic field from equipartition 

The magnetic field Bsa calculated according to Eq. ^ can be 
used to constrain the Doppler factor in an alternative way. The 
magnetic equipartition field Bgq, which minimizes the total en- 
ergy E,oi - (1 + k)Ee + Eg (with relativistic particle energy 
Eg ~ B '-^ and energy of the ma gnetic field Eb ~ B^), is given 
by th e following expression (e.g. lPacholczvld[r970l:lBach et all 
I2005h : 



B,, = [4.5- {I +k)f(a,Va,Vb)LR^) 



2/7 



1.65 ■ 10"' (/(ff, y„, Vi) (1 + z) 5,„ y„, 4 



2/7 



(7) 



with the energy ratio k between electrons and heavy particles, the 
synchrotron luminosity L = An(Pji\ + z) S dv of the source, 
the size of the component R in cm, S ,„ and v,„ in GHz and 
the tabulated function f{a, Va, Vb) with the upper and lower syn- 
chrotron freq uency cutoffs Vg, Vb- Using Vg - 1 0^ Hz, Vb = 5.5 ■ 
lO^Hz (e.g. iBiermann & Strittmatterl 119871: iGhiselhni etall 
IT993I), A: ^ 100 and the same values as for Eq. (|6l), we calculate 
a lower limit for the magnetic field in the range of 1.2-1.3 G, 
which is higher than the magnetic field obtained in the previous 
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section. Eq. (|6]) and O give different dependencies of B and S: 
BsA ~ S, ~ 52/7a+i jYiis yields B^q/BsA = 6%''". Adopting 
the above numbers, we obtain Doppler factors 6eq in the range of 
12-23, in good agreement with both, the previous estimates for 
dvar as well as Svlbi- 



4AA. Equipartition brightness temperature r^'"^ 

Instead of using the IC limit of lO'^ K for the intrinsic brightness 
temperature, a brightness temperature limit based on equipar- 
tition between particle energy and field en ergy could be used 
(IScott & ReadheadLll977tlReadheadill994 : Tg.eg ~ 3 ■ 10" K. 
It is argued that this limit better reflects the stationary state 
of a synchrotron source and for many soii rces is < 10" K 
(e.g. iReadheadL Il994t iGuiiosa & DalvL Il996h . With these num- 
bers we obtain another estimate of the Doppler factor: 6var,eg - 

(1 + z) '"^r7'73 -10" > 8-33, where we again used T"/'' from 

Table|6] Although high, these values still agree with the previous 
estimates and also with Svlbi- 



4.4.5. Inverse Compton Doppler factor 5/c 

Finally we can independently calculate the Doppler factor 
6ic using the upper limits to the soft y-ray flux densities of 
0716+714 obtained from the INTEGRAL observations during 
this experiment and the inverse Compto n argument (see also 
Ostorero et al., 2006; Agudo et al., 2006). Following (see e.g. 
Marscher, 1987; GhiselUni et aL .1993.) the IC Doppler factor 
die is calculated from: 



3/C 



f(a)S„ 



c fl(6-4a) (5-3a) 



1/(4-20-) 



(1+z) 



(8) 



where Vb is the synchrotron high frequency cut-off in GHz, 5,,, 
the flux density in Jy at the synchrotron turnover frequency v,„, 
Sy the observed y-ray flux in Jy at Vy in keV, 9 the source angular 
size in mas and a function o f spectral index f{a) ^ 0.14 - 0.08a 
(e.g. iGhisemni et all Il993h . If we determine the size Ov via the 
causality argument and with the variabi lity time scale fm r, an 
additional 5 comes in Eq. (|2]l (see also lAgudo et all l2006l) and 
we can rewrite Eq. (O: 



(4-2ff)/(10-6a) 



ln{vtlv,n)v^ 



-4a)y(5-3ff) 



1/(10-6q') 



(9) 



whereas now the apparent variability size is calculated from 
= 0.13 tyard^^ (1 + z). Using the same values for v,„, v^, a,/,,,, 
as above and S^ecHz - 5.13 Jy, tyar - 10.7 days (Table |6l), we 
calculate djc using the upper limits to the soft y-ray flux in 
the four bands of INT EGRAL (8, 23, 63 and 141 keV, see also 
lOstorero et al.L |2006|) . The upper limits to the y-ray flux lead 
to lower limits to the IC Doppler-factors, which for the differ- 
ent energy bands a re summarized in Table [T] As outlined by 
lAgudo et al] ( l2006l) . the lower limit of 6ic > 14.1 obtained in 
this way provides a more robust constraint on the Doppler factor 
than Syar and Seq due to the weak dependence of Eq. (|9]l on z and 
di for a,hin - -0.3. We further note that the values of 6jc > 14 
obtained here are in v ery good agreement wi th a Doppler factor 
ofd ~ 17 derived bv lTagliaferri et al.1 (l2003l) . who fit the SEDs 
of 0716+714 by a homogenous, one-zone synchrotron/inverse 
Compton model. 

In Table |7] we summarize the Doppler-factors S calculated 
using the different methods presented above. A comparison of 



Table 7. Summary of the Doppler factor estimates using the dif- 
ferent methods presented in Sect. 14.41 See text for details. 






value 


used values/quantities 


c 

Ovar,IC,ci!! 


> lU—IZ 


ig (laDle[op, ^ gic = i'J 


^varj C,mm 


> 5-10 


Tg (Table|6Jl, rg7c = 10 ^ 




> 8-33 


Tl"' (Table[6ll, Tg"^, = 3 • 10" K 


Seq 


> 12-23 


SsA, S,„ 586G«;=4.0-5.13Jy, 






— (\ 'X Q ^ C\ C\A mac 

^thin — — vJ.j, oyi^^j < u.u+rnas 


SiC.ikeV 


>15.9 


Sikev = 1-37 • 10"' Jy, Q',;„„ = -0.3 






'S86G«z = 5.13 Jy, = 10.7 days 


SlC2ikeV 


> 14.7 


523^^ = 2.55 • 10"^ Jy, a,;,,,, = -0.3 






Sg6GHz = 5.l3}y, tyar = 10.7 days 


SiC.akeV 


>15.3 


5 63fev = 1.23 • 10"^ Jy, ff,;,,,, = -0.3 






'S86GHz = 5.13 Jy, tyar = 10.7 days 


SlCM^keV 


> 14.1 


S 141 kev = 2.42 • 10"' Jy, q',/„„ = -0.3 






'S86G«j = 5.13 Jy, tyar = 10.7 days 



the numbers shows a general good agreement with the Doppler 
factor 6vLBi derived from the superluminal VLBI kinematics. 
Although strong Doppler boosting is required, these values ap- 
pear to be not excessively high, and therefore provide a self- 
consistent explanation of the high apparent variability bright- 
ness temperatures, without a violation of the theoretical limits. 
Since no excessive y-ray emission and therefore no IC catas- 
trophe was recorded by INTEGRAL, and a strong contribution 
of ISS is unlikely, we thus conclude that the violation of the 
theoretical limits inferred from our observations was only ap- 
parent and intrinsic Doppler boosting can naturally explain the 
non-detection of the source at high energies. However, we stress 
that all derived values of S must be regarded as lower limits. 
Consequently, Doppler-factors much higher than 5 > 20 - 30 
can still not be ruled out completely. Contributions from other 
radiation processes, such as partially coherent emission, can not 
be ruled out either 



5. Summary and Conclusions 

We presented the results of a combined variability analysis from 
an intensive multi-frequency campaign on the BL Lac 0716+714 
performed during seven days in November 2003. From seven 
participating radio observatories we obtain a frequency coverage 
of 1.4-666 GHz. Densely sampled IntraDay Variability (IDV) 
light curves obtained at wavelengths of 60, 28, 9, 3, and 1.3 mm 
allowed for the first time a detailed analysis of the source's intra- 
to inter-day variability behavior over the full radio- to short mm- 
band. In this observing campaign 0716+714 was found to be 
in a particular slow mode of variability, when compared to all 
previous IDV observations of this source. While in total in- 
tensity a component of faster variability was observed only at 
/1 60 mm, the source's flux density in the cm- to mm-regime 
was dominated by a nearly monotonic increase on inter-day 
time scales and with variability amplitudes strongly increasing 
towards higher frequencies. Here, our CCF analysis confirms 
that the flux density variations are correlated across the observ- 
ing bands, with variability at shorter wavelengths leading. This 
and the observed frequency dependence, which cannot be ex- 
plained by a model for weak scintillation, strongly suggest that 
the observed inter-day variability has to be considered as source- 
intrinsic rather than being induced by ISS. Only at 60 mm wave- 
length a component of faster variability is seen and implies an 
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unusually high apparent brightness temperature Tb- Hence, ISS 
might be present at this frequency. 

The non-detection of the 'classical', more rapid (type II) IDV 
behavior of 0716+714 in the cm-bands is most likely caused 
by opacity effects and can be related to the overall flaring ac- 
tivity of the source shortly before and during this campaign. 
Sinc e episodic IDV behavior is also observed in other sources 
(e.g. iFuhrmann etall 120021: iKedziora-Chudczed, |2006|) . it ap- 
pears likely that the variability pattern in 'classical' type II IDV 
sources is strongly affected by the evolution of their intrin- 
sic complex structure on time scales of weeks to months. The 
observed complicated variability patterns in total intensity and 
polarization, and their correlations, indicate the existence of a 
multi-component structure with individually varying and polar- 
ized sub-components of different size. 

From daily averages, the spectral evolution of the highly in- 
verted radio-to-sub-mm spectrum of 07 16+7 14 could be studied. 
During the seven observing days the spectrum always peaked 
near v,„ ~ 86 GHz. Significant changes of the peak flux mainly 
occurred during the first 5 days with a continuous rise of the 
peak flux density from 4 to 5 Jy. Together with possibly small 
changes in the daily spectral slopes and the peak frequency v„,, 
the observed variations follow the 'canonical' behavior and indi- 
cate time-variable synchrotron self-absorption and adiabatic ex- 
pansion of a shoc k or a flaring component as described by stan- 
dard models (e.g.l Marscher & Geaflll985h . 

The apparent brightness temperatures Tb obtained from the 
inter-day variations exceed theoretical limits by several orders of 
magnitudes. Although Tb decreases towards the mm-bands, the 
10'^ K IC-limit is always violated. Assuming relativistic boost- 
ing of the radiation, the source must always be strongly Doppler 
boosted. We obtain lower limits to the Doppler factor of the 
source using different methods, including (i) inverse Compton 
(Svarjc) and equipartition (dyar.eg) estimates using the variabil- 
ity brightness temperatures, (ii) an estimate 6^^ using calcula- 
tions of the synchrotron and equipartition magnetic field, and 
(iii) an inverse Compton Doppler factor die using the data from 
the simultaneous INTEGRAL observations. These methods re- 
veal robust and self-consistent lower limits to the Doppler fac- 
tor with 6nir,ic >5, Syar.eq >8 and 5eq > 12. Thcsc limits are 
in good agreement with estimates based on recent kinematical 
VLBI studies of the source and the IC Doppler factor 6jc > 14 
obtained from the upper limits to the high energy emission in the 
3-200 keV bands. The non-detection of the source in the soft y- 
ray bands implies the absence of a simultaneous strong IC catas- 
trophe during the period of our IDV observations. Since a strong 
contribution of interstellar scintillation to the observed inter-day 
variability can be excluded, we conclude that relativistic Doppler 
boosting appears to naturally explain the observed apparent vio- 
lation of the theoretical brightness temperature limits. 
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Appendix A: Statistical variability analysis 

As a criterion for the presence of variability, the hypothesis of 
a constant function (non- variable source) is examined and ax^- 
test of the variability curve S , is done: 



X 



2 

i=i 



( Si- <S > 
[ A5^ 



with the reduced;^^^ 

N 

i 

X7 



A^-l-^l ASi 1 



(A.l) 



(A.2) 



Here, S , denotes the individual flux density measurements at 
time i, < S > their average in time, AS/ the individu al mea- 
surement errors and N the number of measurements (e.g. lTayloil 



119821: iBevington & Robinsoiil Il992h . A source is considered to 
be variable if the -test gives a probability of < 0. 1 % for the as- 
sumption of constant flux density (99.9 % significance level for 
variability). 

To further quantify the observed variability, the modulation 
index m is calculated. This quantity provides a useful measure 
of the strength of the observed variations and is defined as the 
percentage ratio of the rms-fluctuations os of the time series 
and the time averaged mean flux density < S >: 



m[%] = 100 ■ 



<S > 



(A.3) 



In ord er to take the res i dual c alibration errors into account, we 
follow iHeeschen et alj (Il987h and calculate the variability am- 
plitude Y for each individual light curve: 



& Y[%] ^3yfm 



(A.4) 



The F-amplitude can be regarded as being equivalent to a noise- 
bias corrected variability amplitude using the mean modulation 
index mo of the secondary calibrators as a measure of the overall 
calibration accuracy. 

A similar analysis was performed for polarization. The cor- 
responding quantities for Polaris ed intensity S p are 



mp[%] = 100 ■ 



crs, 



<Sp> 



=[%] = 3 ^ 



ml - mj,Q 



For the the polarization angle y we use 
'%[°] = o-^ , 



(A.5) 



(A.6) 



(A.7) 



(A.8) 



A ;t'^-test for S p and x was performed similar to total intensity, 
i.e. using Eq. ( lA.lb and ( IA.2l i. 



Appendix B: Variability time scales 

(i) Structure function analysis 



Following IHeeschen et"aLl (Il987l) and iKraus et all (Il999l) . 
the time series were analyzed by first order structure functions 
(SF). The following procedure was used and applied to each 
data set in order to obtain a variabiU ty time scale at every 
frequency (see also lBeckert et al.L [20021) . For each data set /(?,), 
the discrete structure function 

n 

SF(Tj) = Nij' 2 w(i)w(i + Mm - f{ti + Tj)f (B.l) 

was calculated. Here, Nij is the normalization, f, and f, + Tj are 
the times at which the flux densities / were obtained, and w(i) 
and w(i + j) are weighting functions with w{i)w{i + j) > 0. The 
weighting function accounts for binning of an unevenly sampled 
time series. The binning interval was chosen to be ~ 3-5 times of 
the typical sampling interval of the light curve. In addition, the 
structure functions were also calcul ated using an interpolatio n 
method similar to that described by lEdelson & Krohkl ( Il988h . 
This showed no significant differences in the results compared 
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to the discrete method. Subsequently, two linear fits were per- 
formed to estimate the slope of the SF (SF - a ■ i^) and the 
plateau level (po = 2;«^). This plateau (or saturation) level is 
determined through calculation of the mean value of the SF af- 
ter the first maximum, the related error by the variance of the 
mean value. The formal error of the power law fit to the SF is 
determined by the variance of the fit parameters a and y6. The 
variability time scale is then determined by the crossing point of 
the power law and the fitted saturation level, pQ - a ■ t^, and is 
given by 



tsF 



with its formal error 



AtsF 



a ± Aa/ V3 , 



(B.3) 



In case of no clear saturation of the SF during the whole 
observing period, only a lower limit to the variability time scale 
is given (Table |6] columns 2 and 3). 

(ii) Auto-correlation function analysis 

In the investigation of pulsar scintillation, the usual way is 
to define the "de-correlation" time scale as the time lag during 
which the auto- correlat i on fun ction (ACF) falls to 1/e of its 
maximum (e.g. ICordesL Il986h . The auto-correlation function 
P(t) of a time series S (ti) is defined as tEdelson & Krolik, 198^ 

p(T)^<(S(t)-S(t-T))>, . (B.4) 

The ACF was calculated and tACF was determined by the full- 
width of the ACF at 1/e-height of its maximum. This method 
was applied to our data sets in total intensity and polarization. 
The error estimations were obtained taking the individual 
uncertainties of the calculated auto-correlation functions into 
account. Columns 5 and 6 of Table |6] summarize the resulting 
characteristic time scales and error estimates. 

( Hi) Minimum-maximum method 

The observed total intensity variability patterns are mainly 
characterized by a quasi-monotonic increase of the flux den- 
sity over a significant range of the total observing period. 
Consequently the SF/ACF methods deliver variability time scale 
of 3 - 4 days, which owing to the limited total observing time 
must formally be regarded as lower limits to the true variability 
time scale tyar- In these cases we make use of the following 
definition of the time scale tyar to obtain better estimates of 
the "true" variability time scales directly from those time 
series where only a several day trend was observed (see also 
IMarscher et all 1 19791 and references therein): 



dlnS 1 <S > At 



dt AS 



(B.5) 



' max 



where < S > denotes the mean flux density [Jy] and AS - S,^ 
S ,„i„ the flux density variation in the considered time interval Af . 

We note that this time scale relies on a 'visual' inspection of 
the variability curves, defining by human eye (or mathematically 
via differentiation) the minimum and maximum of the flux den- 
sity and the corresponding time interval Af between these two 
flux density levels. In Table |6] columns 8 & 9 summarize the 
derived values for At and t^,ar of Eq. (|B.5t . 



